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IMPULSIVE MOTION OF SHEAR BUILDINGS INCLUDING 
PLASTICITY AND VISCOUS DAMPING 


Edward Cohen,’ J.M. ASCE, Leon S. Levy? and 
Leonard E. Smollen,* J.M. ASCE 


SYNOPSIS 


The method of normal modes is applied to the analysis of elasto-plastic 
shear buildings subjected to impulsive loads, foundation uplift and non-uniform 
viscous damping. Interaction of foundation rotation and plasticity in the super - 
structure is described, A short table of approximate modes and frequencies 
is included. 


INTRODUCTION 


The use of the normal modes of vibration for the solution of dynamic prob- 
lems involving shear buildings is extended to include the displacement of 
structures subjected to impulsive torsional loads, the development of plastic 
action in the frame and the foundation, ‘‘temporary instability’’ or uplift of the 
foundation and non-uniform viscous damping. The interaction of foundation 
rotation and plasticity in the superstructure is described and a general ap- 
proach to the solution of the problem is outlined. Tables of modes and fre- 
quencies for typical buildings have also been included, 

Shear buildings, as treated in this paper, are defined as idealized multi- 
story structures with the following properties: 


1. The mass of the building is concentrated at the various floor levels. 

2. The motion of the building is due to shear only. The displacements 
due to flexure of the building as a whole are negligible. 

3. The elastic shear in any story is a function only of the elastic de- 
flection between adjacent floors. 


For the purpose of dynamic analysis such idealized buildings are analogous 
to rigid frame structures with stiff floor systems and flexible columns and to 
many other structures of such height to width ratio that overall flexural de- 
formation can be neglected. A typical five story shear building is shown in 
Fig. 1. The equations of motion for a multi-story shear building with impul- 
sive lateral loads will be developed in terms of the normal modes of vibration. 
The use of normal modes is advantageous because the problem is thereby re- 
duced from the solution of the N simultaneous equations that would be required 
by the direct application of Newton's equation of motion to that of N independ- 
ent equations for a building with N degrees of freedom. 

Numerical methods, such as finite difference, trial and error, etc. may be 
used effectively in the solution of the dynamic problems treated in this paper. 
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However, the methods described are also applicable where it is possible to 
obtain analytic solutions. They have been used by the authors for the numeri- 
cal analysis of a three story building subjected to a long duration blast load- 
ing. The computations were done by finite difference theory with the aid of 
only ordinary automatic desk type calculators. However, it is recognized that 
as the number of stories increases the volume of computations also increases 
rapidly and when the number of stories exceeds about five the use of high 
speed digital computers or analog computers becomes necessary. 


A. General Relationships 
The equation of free vibration for the ith floor of a shear building (Fig. 2) 


on a rigid foundation is given by 
Myx, = %) Ky (x - 41) = 1,2,.. (1) 


where 
Mj = The mass concentrated at the ith floor. 
Kj = Spring constant under the ith floor. 
x;(t) = The absolute horizontal displacement of the ith floor at the time t. 


The displacement at any time (t) may be expressed in terms of normalized 
modes ?; and the generalized coordinates, qj(t), as 
N 


By application of the Lagrangian equations, 


4 6T , OU, 


T =the kinetic energy = : 


3 
= the potential energy = — 2 Wj q 
C t P 2 i= J j 


N =the number of stories 

Ww; = the frequence of the jt? mode 

The equations of forced vibration in terms of the generalized coordinates 
become 


N 
ay + = G=1,2,...) (7) 
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. 
(3) 
where: 
N ax; 
= RA P; = the generalized force (4) 
"oe 5 P; = the real force applied to the ith floor 
2 N (5) 
|. 


Thus far it has been assumed that the story masses undergo horizontal 
motion only. This assumption is valid when the foundation material is infinite- 
ly rigid. When the foundation is elastic, the building will also rotate about an 
axis through the centroid of the base and perpendicular to the direction of mo- 
tion as shown in Fig. 3. 
For this case, the following equation of motion may be obtained: 

N 


a, +2 it) = (k= 1,2,...N=1) (8) 


where qx, Qk, and Wj, (i) are the modified generalized coordinates, frequencies, 
and modes of a building rotating elastically about the centroid of its base.‘ 
W,(i) may be expressed as 


N 


where the coefficients C, and Ck; are given by 


a 
A 


%, (i) = Ms the “rocking mode’ of the building 


Ip = the total moment of inertia about the axis of rotation 


The rotation may be expressed as 


(13) 


(14) 


The above equations apply to the case of structures subjected to impulsive 
loads such as may be caused by blast or similar phenomena. The case of 
forced vibrations caused by impressed displacements has been treated 


previously. 


4, M.G. Salvadori. ‘‘Earthquake Stresses in Shear Buildings,’’ Proceedings, 
A.S.C.E., March, 1953, vol. 79,Separate No. 177. 
5. op. cit. 
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| = Cy + Cig H(i; 1,2,..., 
and 
wa (11) 
2k 
N+1 
where 
C 
3, = = Cc, —=— 
Lj y 
Ip 


Ae 


B. Structure Subjected to Impulsive Torque Loading 


When a structure supported on an elastic foundation is subjected to vertical 
impulsive loads having a moment, Mp, (see Fig. 4) about the centroid of the 
base, an additional term must be added to Equation (8). By considering the 
work done by the moment, an additional generalized force due to the 
moment is 


aw i) Cy 
so MrRCk = = MR (1 5) 
Oa, Lj Ip 
and the equation of motion in terms of the generalized coordinates becomes 
oe 2 N 
+Qy = Pj W,,(i) + (16) 


When the moment is caused by vertical loads, and the rotations are small, 
the effect of the equivalent direct force through the centroid and that of the 
couple may be computed independently in the elastic range or as described in 
Section D, when the uplift or plastic action occurs in the foundation. 


C. Plastic Hinges Between Stories in Building on Rigid Foundation 


Plastic hinges are said to develop between any two floors, say the ith and 
the (i+1)t® when the relative semper of the ith and (i+1)*® floor masses 
are such that the shear, K; -Xx i) exceeds the yield point resistance. 
(Fig. 5 shows an idealized i t of Sueer resistance versus displacement). 
From the time a plastic hinge first develops until such time as the relative 
velocities of the two adjacent floors reverse and the story shear returns to 
the elastic range, the method of analysis must take into account the first, and 
all subsequent, plastic hinges as they develop. 

If the problem is solved analytically, it is necessary to revise the equations 
of motion whenever any story exceeds its yield point displacement. This re- 
quirement also holds for a numerical solution. Since the duration of the time 
intervals must be kept to a small fraction *of the period of the given mode 
being computed, the possibility that a new hinge may develop within any one 
step is usually unimportant, the change in equations being made either at the 
beginning or the end of that time interval. 

When considering a shear building in which plastic hinges have developed 
between certain stories, it is clear that the plastic hinges divide the building 
into a number of smaller ‘‘independent’’ structures or substructures, as shown 
in Fig. 6, each with its own modes and frequencies of vibration. Therefore, 
the shear building analysis can be applied to each of the ‘‘substructures’’once 
certain modifications due to the plastic resistance, Sp(i), have been applied 
and the initial conditions established, (Relative motion between the foundation 
and the soil may be treated in exactly the same way.) 

It should be noted that for a building on a rigid foundation the development 
of a plastic hinge between any two floors reduces the problem to the solution 
of two independent ‘‘substructures.”’ Thus, a four story building in which the 


6. The maximum duration of the time interval relative to the modal period is 


a function of the numerical method used and the required accuracy. It may 
also be modified by the variation with time of the loading. 
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third story has become plastic is reduced to two independent 2-story build- 
ings with a consequent simplification of the problem and in the case of numeri- 
cal solutions, a great reduction in the required labor. 

If r stories are isolated by plastic hinges they may be considered to form 
an r-story structure. The equation of motion of the new r-story structure are 
also given by the Lagrange equations, 


d OT OT 


+Q, (k=1,2,... (17) 


where the kinetic energy T, and the potential energy U, are for the new system 
now under consideration. 
The kinetic energy is defined as before, i.e.; 


r 
2 1 
ist (18) 
The potential energy is defined to include the effect of both the internal and 
external work done. This will also include the work done in overcoming the 
plastic resistances, Sp(i), so that once plastic hinges have developed, the po- 
tential energy, U is equal to U, as previously defined, plus the work done by 


the plastic resistances. 


U =U + X(c+r) ~ Sp(e+t) *(c+1) 
or 


The displacements, xj, Can be computed from: 


r 
= 
54 
The initial conditions for the generalized coordinates of any substructure 
isolated by the development of a plastic hinge are determined from the value 
of xP and x? at the instant the plastic hinge is developed. 
At the moment plasticity occurs: 
r 


where the superscript p defines that instant and >; are the modes of the r- 
story ‘‘substructures. 
Multiplying both sides by 95 (i) M, and summing over i, we obtain: 


4 r r 
4 j=l i=l (22) 
However, 
r r 
(i) My = 1 (23) 
j=l i=1 
and 
676-5 


f 
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p 
42, My 


Similarly, 
r 


ap = My 
D. Temporary Instability and Plasticity in the Foundation Soil 


The analysis presented for a rocking building is valid only when the build- 
ing can be considered as rotating about the centroid of the foundation. This is 
approximately true only when the building exerts a compressive force along 
the entire base. However, when the analysis indicates tension on part of the 
foundation the building will pull actually away from the ground, and the center 
of rotation will shift. 

From the soil pressure analysis, a force diagram such as shown in Fig. 7 
may be found. 

The fictitious tensile forces which must be eliminated are equivalent to a 
vertical downward force, T, through the centroid of the base, and a counter- 
clockwise moment, TE, where E is the lever arm of the centroid. An upward 
vertical force, T', and a clockwise moment, T’E, equal in magnitude to the 
quantities to be subtracted are therefore applied to the structure. 

Using a step-by-step finite difference method, a solution by successive 
trials may be obtained as follows: 


1. Assume a value for the vertical force T and the moment TE. 
2. Compute the displacements x; from the equations of Section B by ad- 
justing Equation (16) so that: 


N 
A(i 
+2 a, = 2, (i) + Mp 


where 
(T'E) (i) 


Lj 


dy (27) 


3. Determine the soil pressure diagram based on the values assumed and 
compute T and TE. 

4. If a check is obtained, proceed to the next time interval. If not, repeat 
Steps 1 to 3. 


It may be noted that the proposed method of solution results in a uniformly 
asymptotic approach to the final condition of equilibrium, the computed re- 
sisting moment of the foundation being always greater and the rotation always 
smaller than the actual values. 

Plastic action in the soil may be accounted for in a similar manner by 
equilibrating the difference between the plastic and elastic soil resistance, as 
shown in Fig. 8. 


E. Rotation of Building with Large Plastic Distortions in Superstructure 


Plastic hinges in a building rotating about the center of its base divide it 
into a number of independent substructures, any one of which may be treated 
separately, once the effects of the plastic hinges above and below are included 
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in the analysis. Each plastic hinge is assumed to transmit only a horizontal 
plastic resistance, Spi, an axial force, Tj, and a moment, Mj. Throughout this 
section it will be assumed that the plastic resistance, Spi» is known at any 
given time, t, that it acts in a horizontal direction and may be treated as an 
additional applied force. 

If it is also assumed that the plastic story is axially stiff, i.e., the column 
lengths do not change, both the columns in each story and the floor masses in 
2ach substructure remain parallel and the rotations of all stories are identi- 
cal. The distorted plastic story is then equivalent to a pin connected parallel- 
ogram, or a series of parallelograms when the number of columns exceeds 
two, The adjacent masses move both horizontally and vertically relative to 
each other as the structure rotates. The new unknowns are the horizontal and 
vertical components of the column reactions and the story to story distortion, 
a, as modified by rotation. ~— 

Fig. 9(a) shows the plastic hinge at time, t, where AB is the relative dis- 
placement of the upper substructure at the plastic hinge, anda is the angle of 
distortion of the plastic story. Fig. 9(b) shows the resultant, T, of the final 
axial column forces,and their couple, M, resulting from the rotation. This 
moment is a measure of the contribution of the upper ‘‘substructure’’to the 
total rotational inertia of the structure. The contribution diminishes as the 
distortion, a, increases. The resultant force, T, which is superimposed on 
the pure couple is caused by applied vertical loads and the vertical inertia 
forces. The vertical and horizontal components of the plastic displacements, 
AB, are shown in Fig. 10. 

The following parameters are then obtained for an N-story building with 
n plastic hinges: 


N_ shear displacements x; 


n+1 rotations 6; 


n plastic axial forces Tj 


n plastic displacement angles a ; 


We also have the following independent equations: 


N equations for shear displacement 


n+l equations for rotation 


n equations of rotational compatibility, i.e., 6; = 9), 


n equations of displacement compatibility 


The initial conditions (at rest) are 


x} =0 xf = 0 
=0 5} =0 
a; =0 a; =0 
-O 


T; =0 


For plastic hinge development at any later time, t, suitable initial condi- 
tions may be derived, by the procedure outlined in Section C. 

For a step-by-step solution where the number of plastic hinges is small 
(less than 3 or 4), xj and 6; at the end of any time interval may be found by 
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trial and error by assuming values for a; and T; and computing 6; and xj. The 
equations of shear displacement and displacement compatibility serve as a 
check on the accuracy of the assumed values of aj and Tj. This method is re- 
stricted to a structure with a limited number of hinges because the solution it- 
self has no inherent properties of convergence and depends entirely on the 
computer's judgment and ability to visualize the action of the structure. 


F. Unequal Viscous Damping in Complementary System 


Due to the large difference in magnitude which may exist between the damp- 
ing due to the soil foundation and that due to the shear displacements of the 
structure, it is desirable to account for both these factors, i.e., evaluate the 
actual damping coefficients, My, of the kth mode of the rotating structure in 
terms of the damping coefficients M, of the foundation and/; of the shear 
modes of the structure on a rigid foundation. Both Mj and a may be deter- 
mined experimentally. 

It may be noted that in Equation (1), the equation of motion of any given 
floor mass, viscous damping may be introduced by the addition of linear veloc- 
ity terms, so that Equation (1) becomes: 


N 

x 
where 4;, is the coefficient of the damping between the Ith and ith floor masses 
and where, in general, only the I = (i+1), i, and (i-1) velocities affect the mo- 
tion of the it floor. Therefore 


= Kj 41 - xj) - Ky(xj - + 


The damping coefficients which account for the effect of one story mass upon 
another are difficult to evaluate directly even in the case where rotation is not 
allowed. 

However, if it is assumed that the damping of single degree of freedom 
systems is a linear function of the velocity, a damping term may be added to 
the equation for the generalized coordinates, which becomes: 


N 


The damping coefficient, 4,, may be evaluated in terms of Ay and Hy by con- 
sidering the dissipative function of the damping. 

A dissipative function,’ D, can be defined such that the damping force is the 
negative of the derivative of the dissipative function with respect to the veloc- 
ity. The dissipative function may also be defined by 


1 
D=--W (31) 
2 


where W is the instantaneous work, so that -W is the time rate of energy dis- 
sipation, The dissipative function may be introduced into the Lagrange equa- 
tions, so that: 


7. Karman and Biot, ‘‘Mathematical Methods in Engineering,’’ McGraw-Hill, 
1940, p. 219. 
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d OT OT OV OD (32) 


+ + =0 


dt 04, Oqy, 04, 


Because of the second definition of the dissipative function, the total dis- 
sipative function may be equated to the sum of the soil and shear dissipative 


functions. 

Dy = Dy, + Dyg (33) 
where Dy, the total dissipative function, is: 

N+1 

Dy = Hy (34) 

Dys, the shear dissipative function, is: 
N N N N N Nel 


and Dyg, the soil dissipative function, is: 
N+1 Cc (i) 
2°2 k A 
Dye = “a kat Li 
Substituting in the above equation the values of the dissipative function, and 
differentiating with respect to the kth generalized velocity, we obtain: 


N N 
= Ww 
If we now proceed to the limit where 6, 0 then: 
N N 


Noting that 4; = Hig = + Hi (4-1) and that the 4;;, where I repre- 


sents any story other than the (i+1)™, eal or (i-1), is zero, and knowing H; 
Equation (38) may be solved for the pertinent 4; in terms of uj. These may 
then be substituted in Equation (37) to obtain the required Uy. 


Note: Equation (35) has been used since it is a general mathematical expres- 
sion for a dissipative function in a building of N stories. In physical terms it 
represents the damping action between any two (not necessarily adjacent) 
stories, iandI,. The damping for any floor mass as given by the equation 

for D substituted in the Lagrange equation is equal to OD /Oa, so that the damp- 
ing for the 5t® floor mass, for example, will be given by the following expres- 
sion: 


N N N N 

N N 
- 4%Is 
= 
+ x 
i=] I=) il 5 


5 
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If the effect of damping in the shear modes is small enough to be neglected 
relative to the soil damping, Equation 37 becomes 


My = My 3 (39) 


G. Tables of Modes and Frequencies 


Approximate modes and frequencies of some typical shear buildings having 
2 to 5 stories have been tabulated in Table 1. 


Notation 
M, = mass of the jth story 
Kj = stiffness of the jth story 
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B 
K 
Az 
K) 
| 
+ 
1 


Build 


0.575 5 
1 
=) 0.599 | 0.104 


0.608 


Table I 


M; 
10.75 | 0.615 &0.599 | 0.555 b0.L56 | 0.278 5 1.93 
| 0.559 0.500 L 1.72 
| 0.067 | 0.587 3 1.35 ° 
| 0.496 +0.522 2 | 0.867 
2 | 0.531 | 0.297 1 | 0.298 
| 
core | 
5 |0.40 | 0.637 | -1.03 | 0.835 40.212 | 0.028 | 1.62 
}0.80 | 0.696 0.116 | 0.485 0.115 L 1.29 
3 10.80 | 0.506 | 0.170 0.776 | 0.407 3 0.962 
2 11.00 | 0.322 | 0.552 0.098 40.650 
1 |1.00 |0.150 |0.326 |0.365 |0.551 | 0.660 1 | 0.290 
| 
| 
5 |1.00 |0.599 bo.shé 0322 | 0.165 | 1.92 
L |0.550 0.595 O.L54 1.68 
3 |1.00 | 0-456 10.326 10.550 10.165 | 3 | 1.30 
| 
2 1.00 |0.326 |o.599 |0.17h 10.555 2 | 0.832 
1 |1.00 |0.170 |0.L56 |0.596 |0.331 1 | 0.287 
| | | 
| [1-00 [06658 0.582 $0,233 |1.88 
3 [04577 |0.00 10.582 0.562 3 | 3453 
| 
2 |1.00 lo.ue9 |o.58e 658 2 11.00 
| 
0.228 0.562 |0.658 |o.u29 } 
| 
| | 
| 
| | | 
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1 | 


Table I 


My 
3 | 0.736 0.328 | 3 | 1.80 | 
= 2 }1.00 0.591 | 0.327 L0.736 | | 2 | 1.2L . 
| | 0.327 | | 1 | | 
| | | 
| 
| 2 |1.00 | 0.85b |.0.528 | 2 1.62 
| | 0.528 | 0.851 | | | 1 | 0.616 
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Fig! DISTORTED POSITION OF 
IDEALIZED FIVE STORY BUILDING 
ON A RIGID FOUNDATION 
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Fig. 2 TYPICAL FORCE SYSTEM 
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Xj Lié + Xivs 


O,Axis of Rotation 


Fig. 3 DISTORTED POSITION OF 
IDEALIZED BUILDING ON A_NON- 
RIGID FOUNDATION 
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Moment, M= P-e 


/ (Eccentricity, E | 
Kis M,| / Resultant of Vertical 


Load, 


No. axis of Rotation 


Fig.4 DISTORTED POSITION OF 
IDEALIZED BUILDING ON AN 
ELASTIC FOUNDATION 
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Fig.6 PLASTIC FORCES 
ON "SUBSTRUCTURES 
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Fig. 7 FORCE DIAGRAM ON FOUNDATION FOOTING 
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Fig.8 PLASTIC ACTION IN SOIL 
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Fig.9 REACTIONS AT PLASTIC HINGE 
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Fig.10 VERTICAL AND HORIZONTAL COMPONENTS 
OF PLASTIC DISPLACEMENT 
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592(SA), 593(SA), 594(SA), 595(SA)°, 596(HW), 597(HW), 598(HW)°,599(CP), 600(CP), 601(CP), 
602(CP), 603(CP), 604(EM), 605(EM), 606(EM)°, 607(EM). 


FEBRUARY: 608(WW), 609(WW), 610(WW), 611(WW), 612(WW), 613(WW), 614(WW), 615(WW), 
616(WW), 617(IR), 618(IR), 619(1R), 620(IR), 621(1R)©, 622(1R), 623(1R), 624(HY)©, 625(HY), 
626(HY), 627(HY), 628(HY), 629(H Y), 630(HY), 631(HY), 632(CO), 633(CO). 


MARCH: 634(PO), 635(PO), 636(PO), 637(PO), 638(PO), 639(PO), 640(PO), 641(PO)*, 642(SA), 
643(SA), 644(SA), 645(SA), 646(SA), 647(SA)°, 648(ST), 649(ST), 650(ST), 651(ST), 652(ST), 
653(ST), 654(ST)©, 655(SA), 656(SM)°, 657(SM)°, 658(SM)°. 


APRIL: 659(ST), 660(ST), 661(ST)°, 662(ST), 663(ST), 664(ST)©, 665(HY)°, 666(HY), 667(HY), 
668(HY), 669(HY), 670(EM), 671(EM), 672(EM), 673(EM), 674(EM), 675(EM), 676(EM), 677(EM), 
678(HY). 


c. Discussion of several papers, grouped by Divisions. 
e. Presented at the Atlantic City (N.J.) Convention in June, 1954. 
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